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I. INTRODUCTION 



Nuclear double beta decay is a second order weak interaction process which can proceed 
in two different modes: the two neutrino mode 2v(5(5, with the emission of two neutrinos and 
two electrons, and the neutrinoless mode Qvl3(3 with emission of two electrons only, without 
neutrino emission (see, e.g., Refs. [H, 0, H, 0] ) • Observation of the latter mode that violates 
the lepton number conservation will prove that the neutrino is a massive Majorana particle. 

Theoretical interpretation of electroweak processes in nuclei requires an accurate descrip- 
tion of nuclear many-body wave functions. The proton-neutron quasiparticle random phase 
approximation (QRPA), first considered in Ref. [5j], is one of the most reliable nuclear struc- 
ture methods used for describing the structure of the intermediate nuclear states virtually 
excited in double beta decay. Important ground state correlations are naturally accounted 
for within the QRPA as well (see, e.g., Refs. 0, i, 9). The method has been shown to 
be capable of successful describing double beta decay provided the particle-particle (p — p) 
residual interaction is included into consideration, along with the usual particle-hole {p — h) 
one [3]. On the other hand, the calculated nuclear matrix elements for double beta decay, 
both the two-neutrino and neutrinoless modes, have been proved to sensitively depend on 
the parameter Qpp that renormalizes the G-matrix strength in the p — p channel [l, 1^. 

As the majority of /5/3-decaying nuclei are nearly spherical, the spherical symmetry has 
been assumed in all QRPA calculations of the Oz//3/3-decay matrix elements M^'^ up to now. 
Nowadays, there is a growing interest in the scientific community in double beta decay of 
^^°Nd as a very promising candidate for the next generation of experimental searches for the 
neutrinoless double beta decay (SNO-f- [sl and SuperNEMO [^). The nucleus is well known 
to have one of the largest values of the lepton phase space factor (about 33 times larger than 
that for ''^Ge, see, e.g. Ref. but it is at the same time strongly deformed that provides 
a great obstacle for reliable theoretical analysis of the corresponding matrix element M°^. 
The effect of deformation on the Oz//3/3-decay of ^^°Nd to the ground and excited states of 
^^°Sm was first studied in Ref. [lO| within the pseudo-SU(3) model that is a tractable shell 
model for deformed nuclei, in which, however, a severely truncated single-particle basis is 
used. The matrix element M^" = —1.57 was found in Ref. [lo| for the ground state to 
ground state (g.s.-to-g.s.) transition ^^°Nd— >-^^°Sm, substantially smaller than the recent 
result M^'^ = 4.74 of Ref. [6|. In the latter publication, however, the nuclei were considered 
as spherical. Recent result M^'^ = —1.61 obtained in Ref. ll| within the Projected Hartree- 
Fock-Bogoliubov approach is in a good agreement with the previous one of Ref. El. Also 
the double beta decay ^^°Gd— >^^°Dy is of interest although the lepton phase space factor is 
about 6 times smaller than in the case of ^^^Nd^^^'^Sm. 

An extension of the pnQRPA method to accommodate the effect of nuclear deformation 
was first done in Ref. [l2|, where single particle (s.p.) basis was generated in a Nilsson 
potential. Further developments by including Woods-Saxon type potentials [l^], residual 
interactions in the particle-particle channel [l4l| , self consistent deformed Hartree-Fock mean 
fields with consistent residual interactions |l5l| have followed the article [l^ . 

Recently, the Gamow- Teller strength distributions and the 2z//3/5-decay matrix elements 
have been calculated for deformed nuclei within the QRPA by making use of a phenomeno- 



logical deformed Woods-Saxon potential and the schematic separable forces [16|, ll7|, Il8| . It 
has been found that differences in deformation between initial and final nuclei can have a 
pronounced effect on the 2z//3/3-decay half-lives. However, using the schematic forces for cal- 
culating the amplitudes of the Oi//3/?-decay would immediately raise the problem of how to 
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fix the numerous strength parameters of the forces in different J'^ partial channels. There- 
fore, using realistic interaction with a minimal number of renormalization parameters is 
obviously preferable. In the present paper the approach used in Refs. [13, lH| is extended to 



accommodate realistic effective interaction based on the Brueckner G-matrix derived from 
the nucleon-nucleon Bonn CD force. 

The 2i//3/3-decay half-lives have been already measured for a dozen of nuclei and the corre- 
sponding nuclear matrix element M^^^ have been extracted [l9|. Theoretical interpretation 
of these matrix elements provides a cross check of reliability of the calculated nuclear wave 
functions. In recent most detailed spherical QRPA calculations of the Oi//5/3-decay matrix 
elements Ref . [6[ , the renormalization parameter Qpp of the particle-particle residual inter- 
action is fixed in such a way that the experimental half-lives of 2z//5/3-decay are correctly 
reproduced. Thus, before calculating the matrix elements M^^ using the QRPA extended 
to deformed nuclei, the corresponding calculation of M"^" has to be done as a test for the 
modeled nuclear wave functions. 



II. THE QRPA IN DEFORMED NUCLEI 

The inverse half-life of the 2i//3/?-decay can be expressed as a product of an accurately 
known phase-space factor G^^ and the second order Gamow- Teller matrix element M'qj^ for 
the g.s.-to-g.s. transition 

Contribution from the two successive Fermi transitions to the amplitude of the 2z//3/5-decay 
can safely be neglected as it arises from isospin mixing effects (see, e.g., Ref. jl]). The double 
Gamow- Teller matrix element in Eq. ([1]) can be written in the following form: 

,,._j,(ofl^1WMm), (2) 

where the index i{f) refers to the initial (final) nuclei, the sum runs over all |m) = ll"*") states 
of the intermediate odd-odd nucleus, /?~ = ^oJa the Gamow- Teller transition operator, 

LOm = Em - (-Eo, + -EoJ/2 = {u)m{i) + U;m(/))/2 with U)m(i) = - Eo^ {i^m{f) = Em - Eq^.) 

representing the excitation energy of the m'th state relative to the g.s. of the initial (final) 
nucleus. 

To take into account the effect of deformation, wave functions 11"*") of the intermediate 
states in the laboratory frame which have a projection M of the total angular momentum 
onto the axis z can be represented in terms of wave functions in the intrinsic frame: 



\lM{K),m) = ,l^,[Vl,A<P,0,^P)Ql^+{-lY^^Vl,_^^^^ {K = ±1), 

\lM{K),m) = J^Vl,MO,ml,KKs.) iK = 0). (3) 



Here, the correlated QRPA ground state in the intrinsic frame is denoted as 10^^ ), the 

intrinsic excitations are generated by the QRPA phonon creation operator Qm^^ ^ 
the projection of the total angular momentum onto the nuclear symmetry axis (the only 
projection which is conserved in strongly deformed nuclei). 
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These adiabatic Bohr-Mottelson-type wave functions provide an approximation which is 
vahd for large deformations. Thus, the most interesting nuclei in question, ^^°Nd and ^^°Sm, 
which indeed are strongly deformed, are well justified to be treated in such an approximation. 
The adiabatic, or the strong coupling (see, e.g., the monograph [sO]), approach fails, however, 
for small deformations since the Coriolis force gets large and mixes states with different K. 
Nuclei ''^Ge and ™Se have rather small deformations and the so-called weak coupling, or no 



alignment, limit [20| seems to be more suitable. In this limit the Coriolis force becomes so 
strong that the angular momenta of the valence nucleons get completely decoupled from the 
orientation of the core. Such a case would deserve a detailed study that is postponed to a 
future publication and the adiabatic approach to description of excited states of deformed 
nuclei is adopted in the present first application of the QRPA with a realistic residual 
interaction. Nevertheless, one might already anticipate without calculations that in the 
weak coupling limit the calculated observables should reveal smaller deviations from the 
ones obtained in the spherical limit than those calculated in the strong coupling limit of 
the present work. In this connection it is worth to note that spherical QRPA results can 
exactly be reproduced in the present calculation by letting deformation vanish, in spite of 
the formal inapplicability of the strong-coupling ansatz for the wave function in this limit. 
The QRPA phonon creation operator acting on the ground-state wave function is given 

as: 

Qm.K = y^.^pn.K^pn.K ~ ^pn.K^pn,K- (4) 
pn 

Here, J^^n^K — (^\(^n and Apn,K = a-pO^n are the two-quasiparticle creation and annihilation 
operators, respectively, with the bar denoting the time-reversal operation. The quasiparticle 
pairs pn are defined by the selection rules f2p — f2„ = i^' and 7rp7r„ = 1, where vr^ is the single- 
particle (s.p.) parity and is the projection of the total s.p. angular momentum on the 
nuclear symmetry axis (r = p, n). The s.p. states |p) and |n) of protons and neutrons are 
calculated by solvin g th e Schrodinger equation with the deformed axially symmetric Woods- 
Saxon potential jli, 2^. In the cylindrical coordinates the deformed Woods-Saxon s.p. wave 



functions \t^t) with $1,- > are decomposed over the deformed harmonic oscillator s.p. 
wave functions (with the principal quantum numbers (A^n^A)) and the spin wave functions 
|E = ±1): 

\r^r) = E ^^n.s|A^ri,A. = f^r-S)|S), (5) 

where A^ = n_L + (n_L = Inp + |A|), and Up are the number of nodes of the basis 
functions in the z- and p-directions, respectively; A = f2 — E and E are the projections 
of the orbital and spin angular momentum onto the symmetry axis z. For the s.p. states 
with the negative projection = — I^^t-I, that are degenerate in energy with = \^t\, the 



time-reversed version of Eq. ([5]) is used as a definition (see also Ref. [17|). The states (t, f) 
comprise the whole single-particle model space. 

The deformed harmonic oscillator wave functions \NnzA) can be further decomposed 
over the spherical harmonic oscillator ones In^/A) by calculating the corresponding spatial 
overlap integrals ^ = (nj./A|A^n^A) (n^ is the radial quantum number, I and A are the 
orbital angular momentum and its projection onto 2;-axes, respectively), see Appendix for 
the details. Thereby, the wave function ([5]) can be reexpressed as 

\rn,) = Y.B:,\v^r), (6) 
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where Ir/fi,) = ECf^^^^ i ^| 



nJA = — S)|E) is the spherical harmonic oscillator wave 



function in the j-coupled scheme (77 = (rirlj)), and BJ^ = J^C] 
with C'!^ ^ 1 ^ being the Clebsch-Gordan coefficient. 



The QRPA equations: 



A{K) B{K) \ 
-B{K) -A{K)) 



Xm 
K 



K 



Xm 
K 
Vm 



(7) 



with realistic residual interaction are solved to get the forward XJ^, backward F-^ ampli- 
tudes and the excitation energies cu'^' and co^^ of the m-th {K = 0, ±1) state in the 
intermediate nucleus. The matrix A and B are defined by 



Apn^p'n' 



^pn,p'n' 



{Ep + En) + gppiUpUnUp'Un' + VpVnVpiVn')Vpnp'n' 
- gphiUpVnUpiVni + t^pM„fp'M„' ) ^pn'p'n 
dppiUpUnVpiVn' + VpVnUpiUn')Vpf^p/^/ 
QphiUpDnVpiVn' ~\~ 'VpUnUp'Vn')Vpn'p'n 



(8) 



where Ep + En are the two-quasiparticle excitation energies, l^„y„' and Vpn^p,-, are the p — h 
and p — p matrix elements of the residual nucleon-nucleon interaction V, respectively, Ur 
and Vr are the coefficients of the Bogoliubov transformation. The amplitudes of /3~ and 
f3~^ transitions from the O"*" g.s. of initial and final nuclei to a one-phonon state in the 
intermediate nucleus are given in the intrinsic system by: 



UpVnX'!^^^ + "VpUnY^j^ 



pn 



(K+,m|/?+|o+,.) = Y.(pWK\n) lupVnY;;^^^ + vpUnX;;n,K 



(9) 



pn 



The matrix element Mq^ of Eq. ([2]) is given within the QRPA in the intrinsic system by 
the following expression: 



E E 

K=0,±1 mi-nif 



{Om],\K+,mj){K+,mj\K+,m,){K+,m,\P],\Ot) 



(10) 



Hmf 



mimf 



Along with the usual approximation of the energy denominator in Eq. ( ITOl) as uk. 
{^K,mf+^K,mi)/'^ (sec, e.g., Refs. [l7|,[l8|; we will later refer to this case as "case 11"), we also 
use in this work another prescription in which the whole calculated QRPA energy spectrum 
is shifted in such a way as to have the ffist calculated 1^ state exactly at the corresponding 
experimental energy (case I). In this case the energy denominator in Eq. (fTOj) acquires the 
form ujK,mimf = {ooK,mf " i^K,if + ^K,m, " ujk,u)/'2 + co'i+ , with 00-^+ being the experimental 
excitation energy of the ffist 1+ state measured from the mean g.s. energy [Eq. + £'o^)/2. 
All the calculated strength functions in this work are represented according to the case I, as 
well. 

The two sets of intermediate nuclear states generated from the initial and final g.s. do 
not come out identical within the QRPA. Therefore, the overlap factor of these states is 
introduced in Eq. (fTOj) 22, H] as follows: 



{K+, mf\K+, m,) = [X^^'kKk ' YiJ^yi"^] ^i.k {BCSf\BCS,). 
kif 



(11) 
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The factor TZiji^, which includes the overlaps of single particle wave functions of the initial 
and final nuclei is given by: 

nw = (j9p,bV)(«?«P'^ + v^^vlP){np^\n'pr.,){u^^ui^,^ + v^^vlP), (12) 

and the last term {BCSf\BCSi) in Eq. (ITTl) corresponds to the overlap factor of the initial 
and final BCS vacua in the form given in Ref. [l7| . 

As a residual two-body interaction we use the nuclear Brueckner G-matrix, that is a 
solution of the Bethe-Goldstone equation, derived from the Bonn-CD one boson exchange 
potential, as used also in the spherical calculations of Ref. The G-matrix elements are 
originally calculated with respect to a spherical harmonic oscillator s.p. basis. By using the 
decomposition of the deformed s.p. wave function in Eq. ([6]), the two-body deformed wave 
function can be represented as: 

\pn) = i^/,f„,J^pr/„,Jir), (13) 

where \vpVn, JK) = C^ii,„nJ^P^^p) l^n^n), and = i?,^,i?,^(-l)^"-^"Cii,„_o„ is 

defined for the sake of simplicity ((—1)-'" ^" is the phase arising from the time-reversed states 
\n)). The particle-particle Vp^^ p,-, and particle-hole Vpn\ p'n interaction matrix elements in 
the representation for the QRPA matrices A, ([7]) in the deformed Woods-Saxon single- 
particle basis can then be given in terms of the spherical G-matrix elements as follows: 

Vpn, p'n' =~ ^ J^pil^nrir,^p'^p'n''n„t^i''lpVnVp'Vn', J), (14) 

J VpVn Vp'Vn' 

Vpn',p'n = ^ ^pr,pn"v„'^p'Vp'mn'^iVpVn'Vp'Vn,J), (15) 

J VpVti Vp'Vn' 

where Kp^, = Qp + Qn' = ^p' + ^n- The matrix elements of ax in Eq. can be written as 
{p\aK\n) = j:^p,nnPp^pnvJVp\W\\Vn)/VS- 



III. RESULTS 

The Gamow- Teller strength distributions and the 2i//?/9-decay amplitudes for the g.s.-to- 
g.s. transitions are calculated for the nuclear systems with A = 76 (^^Ge — >™Se) and A 
= 150 (^^°Nd — i>^^°Sm). The single-particle Schrodinger equation with the Hamiltonian of 



a deformed Woods-Saxon mean field [24, |25j is solved on the basis of a axially-deformed 



harmonic oscillator, see Eq. ([5]). The s.p. basis corresponding in the spherical limit to full 
{2-4:)huj major oscillator shells for the nuclei with A = 76 and (4-6)^0; for A = 150 are 
used. Decomposition over the states within seven major spherical harmonic oscillator shells 
is done in Eq. (fT3l) . Only quadrupole deformation is taken into account in the calculations. 
The deformation parameter P2 is obtained as P2 = \f^'§^ (^c is the charge rms radius) by 
using the empirical intrinsic quadrupole moments Qp which are derived from the laboratory 
quadrupole moments measured by the Coulomb excitation reorientation technique [2l|. The 
corresponding values of P2 are listed in Table [H The spherical limit, i.e. P2 = 0, is considered 
in all cases as well. Because of rather large experimental errors (32 = 0.37 ± 0.09 and 



j32 = 0.23 ± 0.03 [2l| for ^^°Nd and ^^°Sm, respectively, we also adopt for these nuclei the 
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28 which seem to fit better the rotational bands in 



respective calculated values from Ref. 
these nuclei. 

First, the BCS equations are solved self-consistently to obtain the occupation amplitudes 
Ur and Vr, gap parameter and the chemical potentials Xp and A„ 2^. The renormalizing 
strengths ^f^^j^ and (7^^^^ of the proton and neutron pairing interactions are determined to 
reproduce the experimental pairing energies through a symmetric five term formula (27| . 

For calculating the QRPA energies and wave functions one has to fix the particle-hole Qph 
and particle-particle Qpp renormalization factors of the residual interaction in Eq. ([8]). An 
appropriate value of gph can be determined by reproducing the experimental position of the 
Gamow- Teller giant resonance (GTR) in the intermediate nucleus, whereas the parameter 
Qpp can be determined from fitting the experimental value M^^"*^^^. The experimental energy 
position of the GTR relative to the energy of the first excited l"*" state in ^^Ge can be repro- 
duced with gph = 1.15 (the prescription of case I is used, see explanations after Eq. flTUl) ). 
Since there is no experimental information on the GTR energy for ^^''Nd, we use for this 
nucleus the same gph = 1.15. Two sets of parameters gpp obtained from fitting M^rf^^^ in 
two calculations of Mq^j, (case I and case II) are listed in Table [B The difference between 
gpp fitted in case I and case II is usually quite small. To compare with the previous QRPA 
results, we have performed calculations also with the separable p — h and p—p interactions in 



a way similar to what was done in Refs. [17|, Il8| . The coefficient in the A-dependence of the 



p — h strength parameter x = 3.73/ A^-'^ MeV is fitted to reproduce the GTR energy in ''^Ge, 
and is used then in the calculations for ^^°Nd (the form of the A-dependence is taken from 
Ref. The corresponding fitted values of the strength parameters k of the separable 

p — p interaction are also listed in Table [B All the calculations in this work are done with 
the mean field spin-orbit coupling constant increased by factor 1.2 as compared with the 
one used in Refs. 17, 1^, that gives a better correspondence with the parametrization of 
the spherical Woods-Saxon mean field used in Ref. m ■ Therefore, we get slightly different 



fitted values of x and k as compared with Ref. jl7l . [l 



TABLE I: The values of the deformation parameter /32 for initial (final) nuclei adopted in the 
calculations along with the fitted values of the p — p strength parameters gpp (for the realistic 
Bonn-CD force) and k (for a phenomenological separable force) for two ways of calculations of 
Mo^rp, (I) and (II). The p — h strength parameters gph = 1-15 and x = 3.73/j4'''^ MeV are fixed as 
explained in the text. In the last column the calculated (with the Bonn-CD force) values of the 
Ikeda sum rule (in % of 3(A^ ~ ^)) are given for the initial (final) nucleus. 



nucleus 


(32 


9pp (I) 


9pp (II) 


K (I), MeV 


K (II), MeV 


ISR(%) 


76Ge C^Se) 


0.0 (0.0) 


0.94 


0.91 


0.087 


0.083 


96.8 (98.2) 




0.10 (0.16) [21] 


0.99 


0.97 


0.091 


0.088 


96.1 (96.8) 


iso^d (isogm) 


0.0 (0.0) 


1.11 


1.11 


0.051 


0.050 


94.8 (95.9) 




0.37 (0.23) [21] 


0.78 


0.65 


0.033 


0.005 


94.1 (95.8) 




0.24 (0.21) [28] 


1.35 


1.32 


0.053 


0.052 


95.4 (96.2) 



The calculation of single (3~ and (3~^ decay branches for parent and daughter nuclei is the 
starting point for the calculation of the 2z//3/3-decay amplitudes. The calculated Gamow- 
Teller strength distributions for all the nuclei in question are shown in Figs. 1,2 as functions 
of the excitation energy in the intermediate (for the 2z//3/3 decay) nuclei. The representation 
is according to the convention of the case I, i.e. the entire calculated GT spectrum is shifted 
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in order to fit the experimental energy of the first l"*" state. In order to facihtate comparison 
among various calculations the Gamow- Teller distributions are smoothed with a Gaussian 
of width 1 MeV, so the original discrete spectrum of B{GT) values is transformed into 
a continuous one of the strength function S{GT). The relevant strength distributions for 
the double beta decay are the S{GT~) for the parent nucleus (two upper panels) and the 
S{GT~^) distribution for the daughter nucleus (two lower panels). In the left panels labeled 
as "Realistic" the results obtained by using a realistic nucleon-nucleon interaction (Bonn 
CD force) are shown whereas in the right ones labeled as "Separable" the results obtained 
by using a separable interaction of Ref . [3] • In both types of calculations the corresponding 
strengths of the p — h interaction Qph and x ^'^'^ fixed by fitting the experimental energy of 
the Gamow- Teller resonance in ^^Ge as described in the preceding paragraph. The thick 
solid and the dashed lines represent the results obtained in the spherical limit = and 
with the realistic deformation, respectively. The thin solid line in Fi g. 1 represents the 
Gaussian-smeared experimental B(GT) values for ^^Ge taken from Ref. 29l . 

In the case of I3~ distribution, one observes that the position of the Gamow- Teller reso- 
nance is not sensitive to the effect of deformation. In the case of (3^ distribution, the effect 
of deformation is more apparent in the strength distribution than that in the case of the j3~ 
one. Comparing the results obtained with the realistic and schematic residual interaction 
one can see some marked differences in the and the low-energy part of the j3~ strength 
distributions. The Ikeda sum rule is underestimated by a small amount about (3-5)% in 
the calculations (see last column of Table I; to get 100 % one would need to have the whole 
s.p. basis, as it is done, for instance, within the continuum-QRPA, see e.g. Ref. [13]). It 
also means that the chosen s.p. model spaces (see the beginning of this section) are large 
enough. 

Fig. 3 illustrates the evolution of the strength functions S{GT^) in '^^Ge and ^^'^Nd with 
respect to increase of Qpp in the spherical limit {(32 = 0). One can see in the figure that all 
the GT peaks get shifted to smaller energies as gpp increases. In addition, the low-lying GT 
states become more collective and, correspondingly, the low-energy part of the GT strength 
gets markedly larger at the expense of a decrease in the GTR strength. It can also be seen 
that the QRPA calculations with the realistic Qpp as listed in Table I are still quite far from 
the collapse of the QRPA. One may argue that the results show too enhanced low-energy 
part of the GT spectrum as compared with the experiment. It might have to do with neglect 
of quenching in the present calculation. For qa < 1.25 the corresponding experimental value 
of M^rf^^^ gets larger and a smaller Qpp would be needed to fit it. This effect deserves a 
separate detailed study, similar to the one performed in Ref. j3l|, that is out of scope of the 
present work. 

As shown in Ref. [l7, 18], deformation introduces a mechanism of suppression of the 
matrix element that works even for the same initial and final deformations. The reduction 
gets even stronger when initial and final deformations differ from each other, and is mainly 
due to a decrease of the BCS overlap (ITTl) . The values of the overlap calculated with the 
realistic pairin g in teraction are very close to those obtained by using the constant pairing 
gap as in Refs. 17|, ll8|] that is illustrated in Fig. 4 for the transition ''^Ge— >-^^Se. In addition, 
there is the well-known reduction of the with increasing renormalization parameter Qpp 
of the p — p interaction 0, 0]- It is interesting to study all these suppression mechanisms 
within the QRPA with the realistic residual interaction that is main thrust of this work. 

The calculated matrix elements are shown in Figs. 5,6 as functions of the parameter 
Qpp and for different deformations of initial and final nuclei. The parameters j32, Qph and 
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X are used the same as in the study of the Gamow- Teller distributions. Again, in the left 
panels labeled as "Realistic" the results of present work using the realistic nucleon-nucleon 
interaction (Bonn CD force) are shown whereas in the right ones labeled as " Separable" the 
results using separable interaction are presented. The upper panels represent the results 
obtained with the shifted calculated QRPA spectrum (case I, see explanation after Eq. (ITUil ) 
and in the lower ones the results obtained with the usual unshifted QRPA spectrum (case II) 
are shown. The solid lines in Figs. 5,6 represent the matrix elements M'^rp calculated in the 
spherical limit while the dashed ones (and dot-dashed in Fig. 6) represent the M'^j- calculated 
for realistic deformations. The dotted horizontal line corresponds to the corresponding 
experimental values Mq^^^^^ obtained in Ref. [ii by using unquenched value of the axial- 
vector coupling constant qa = 1.25. The points A and B (and C for ^^^Nd^^^'^Sm) in each 
panel specify the values of the p — p interaction for which the value of M^^"^^^ is reproduced 
for spherical and deformed cases, respectively (the corresponding values of g^p and n are 
listed in Table I). The calculated decreases faster in the spherical case than in the 

deformed one as Qpp increases, and the fitted value of Qpp is usually larger in the case of 
deformation than that for the spherical limit. Another interesting feature is that for non- 
zero deformation the calculated M^^^ smaller at q^n = than those in the spherical limit, 



in agreement with the previous results of Refs. |17l. Il8[|. This suppression already holds when 



initial and final nuclei are equally deformed and gets stronger with increasing difference in 
deformation between initial and final nuclei (the BCS overlap plays an important role in 



the reduction of M^.^^, as discussed in Refs. [17|, Il8[). It is also worth to mention that 
the calculated M^^ for i^°Nd^^^°Sm transitionjTig. 6) comes out rather small when the 
experimental deformation parameters from Ref. 2l| are used, while by using the deformation 
parameters of Ref. 28] the corresponding M^rf^^^ can be fitted for a reasonable value of 
Qpp. Also, the value of M^j, calculated with the shifted QRPA spectrum, though a little 
bit larger at gpp -C 1 than the unshifted ones, fit the corresponding experimental values at 
almost the same gpp (see Table I). 



IV. CONCLUSIONS 

In the present work the two- neutrino double beta decay (g.s.-to-g.s. transitions) and rel- 
evant Gamow- Teller strength distributions are calculated within the QRPA for the nuclear 
systems ^^Ge— >^^Se and ^^°Nd— ^^^''Sm with taking into account effects of nuclear deforma- 
tion. For the first time a realistic residual two-body interaction based on the Brueckner G- 
matrix (for Bonn CD force) is implemented in deformed calculations. The G-matrix elements 
in the deformed Wood-Saxon basis are calculated by expanding the deformed single-particle 
wave functions over the spherical harmonic oscillator ones. The effect of deformation on the 
observables is studied within the framework. The suppression of the calculated M^^^, due to 
both non-zero deformation and the interaction in the particle-particle channel, is observed 
in accordance with previous calculations with separable forces (itI . 18| . The present work 



reahzes the first important step towards the QRPA calculations of the neutrinoless double 
beta decay of deformed nuclei like ^^°Nd with a realistic nucleon-nucleon interaction. 
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APPENDIX A: CALCULATION OF THE OVERLAP INTEGRALS 



In this Appendix we describe how the deformed harmonic oscillator wave functions can 
be decomposed over the spherical ones by calculating the corresponding spatial overlap 
integrals. 

The normalized wave equation of the three-dimensional axially-deformed harmonic oscil- 
lator in cylindrical coordinates (p, 0) are given by a product of three functions: 

|A^^.A)=<'(p)V'n.(^) (Al) 



27r 



where A the projection of the orbital angular momentum on z (symmetry) axis, N is the 
principal quantum number which is defined as = n^; -|- 2np + |A]. The radial function -01^1 
is usually written in terms of a dimensionless coordinate 77 as 

<l(p) = ^^]il^^e-iLlfl(r;) (A2) 



with ?7 = where h±_ = \J -^^j^ oscillator length for the motion perpendicular to the 

z-axis, C]^^ = +|A|)!b^ ^ ^ normahzation factor and L\^^{r]) are the associated Laguerre 

polynomials. The z-dependent function ip^^ is similarly written in terms of a dimensionless 
variable ^ as 

^r^M-Cr^^e-^H^M (A3) 
with ^ — here is oscillator length in the direction of the z-axis, = (y^2"'^n2! bz)~^ 
is a normalization factor and HnXV) ^'^^ the Hermite polynomials. 

The normalized wave functions of the three-dimensional isotropic harmonic oscillator in 
spherical polar coordinates (r, 6*, 0) are given: 

\NlK)^i;r,^i{r)Yi^{e,4>), (A4) 

where rir is the radial quantum number, Yi\ is the spherical harmonic, / is the orbital angular 
momentum. The radial part il^nri is written in terms of a dimensionless coordinate v 

^l(r)^Cr.y/'e-^df^\u) (A5) 

with i/ = where 60 is the spherical oscillator length, Cnj — (7 — _^JYu,u3 ) is a normal- 
ization factor and Lur ^ (i^^) are the associated Laguerre polynomials. 
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The wave functions |A^n^A) of the deformed harmonic oscillator can be decomposed over 
the spherical ones \NIA) as 

\NnA) = E KLaINIA) (A6) 

where ^IjvLa ~ {NlA\NnzA) is the spatial overlap integral which can be numerically calcu- 
lated in the spherical coordinate system as follows: 

AnLa = ^ rf/V|fJ(rsin^?)^„^(rcos^)F,l(^,0 = O)sin^^(i^?)^,,K^)^'c?r. 



vo 
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FIG. 1: Gamow- Teller strength distributions S{GT^ ) m 
of the excitation energy E^x in the intermediate (for ^^Ge^''^Se decay) nucleus ^^As. The QRPA 
calculation results obtained with realistic and separable forces are shown in panels (a) and (b), 
respectively. Upper panels: the results corresponding to the spherical (/32 = 0.0) and deformed 
{P2 = 0.1) ground state of ^^Ge are represented by the thick solid and dashed lines, respectively. 
Experimental data (thin solid line) are from [29t |. Lower panels: The results corresponding to the 
spherical {(32 = 0.0) and deformed {(32 = 0.16) ground state of ''^Se are represented by the thick 
solid and dashed lines, respectively. 



13 



Realistic 



Separable 




11 — cal. , ^=0.0 
- - cal. , p2=0-37 __ 
... cal., p2=0.24 — 



cal., ^2=0-0 
cal., ^=0.23 
cal., ^2=0-21 



10 15 20 

E [MeV] 




25 







10 



15 



20 



25 







30 



E [MeV] 

ex. 



FIG. 2: The same as in Fig. 1, but for ^^^Nd and ^^^Sm. 
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FIG. 3: Gamow- Teller strength distributions S{GT~) in Ge and Nd (spherical limit) for 
different particle-particle interaction strengths gpp as functions of the QRPA energy Eqrpa- The 
solid lines (the same as in Fig. 1,2, upper left panels) represent S{GT~) corresponding to the fitted 
gpp (Table I). The dashed and the dotted lines correspond to smaller gpp as indicated in the Figure. 



15 



Realistic 



Separable 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




_ p2(Ge)=-0.3 
p2(Ge)=-0.2 7 
p2(Ge)=-0.1 
p2(Ge)=0.0 
p2(Ge)=0.1 
p2(Ge)=0.2 
p2(Ge)=0.3 




-0.3-0.2-0.1 0.1 0.2 0.3 
P2(Se) 



-0.3-0.2-0.1 0.1 0.2 0.3 

P^CSe) 



FIG. 4: The dependence of the BCS overlap factor in Eq. (jlip on the difference in deformations 
of ^^Ge and '^^Se. The calculation results obtained with realistic and separable forces are shown in 
panels (a) and (b), respectively. 
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FIG. 5: 2z^/3/3-decay matrix element for Ge^ Se decay as function of particle-particle interaction 
strength Qpp of the realistic forces (left panels), and of n of the separable forces (right panels). 
The solid line ("sph.") corresponds to spherical shape of the initial and final nuclei. The dashed 
one ("def.") is associated with the deformation parameters taken from Ref. [2l| (/32('^^Ge)= 0.1, 
/32(^^Se)= 0.16). The doted horizontal line corresponds to experimental M^^*^^^ obtained in 
Ref. [19] using qa = 1.25. The points A and B in each panel specify the values of the p — p 
interaction for which the value of M^^*^^^ is fitted. The upper (case I) and lower (case II) panels 
correspond to the calculations with the shifted and unshifted QRPA spectrum, respectively, as 
explained in the text. 
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FIG. 6: The same as in Fig. 5, but for ^^'^Nd- 
parameters: from Ref. (/32(^^°Nd)= 0.37, /32C^"Sm 
(/32(^^°Nd)= 0.24, /?2(^^°Sm)= 0.21, "def. (2)"). 



150 decay with 2 different sets of deformation 
0.23, "def. (1)" ) and from Ref. |2£l 



18 



